Influence of oxygen ordering kinetics on Raman and optical response in YBa 2 Cu 3 06.4 
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Kinetics of the optical and Raman response in YBa2Cu306.4 were studied during room temperature 
annealing following heat treatment. The superconducting T c , dc resistivity, and low-energy optical 
conductivity recover slowly, implying a long relaxation time for the carrier density. Short relaxation 
times are observed for the B\ g Raman scattering - magnetic, continuum, and phonon - and the 
charge transfer band. Monte Carlo simulations suggest that these two relaxation rates are related to 
two length scales corresponding to local oxygen ordering (fast) and long chain and twin formation 
(slow) . 
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Room temperature (RT) annealing has been observed 
in oxygen deficient YBa2Cu306+x single crystals rapidly 
quenched from elevated temperatures (T w 420 K) JjJ. It 
is generally believed that the disorder produced by heat- 
ing is frozen by the fast quench, and that RT anneal- 
ing results from re-ordering of oxygen in the CuO-chain 
layer via diffusion. The oxygen re-ordering lowers the 
average valence of the chain Cu and transfers holes to 
the Cu02-planes. Thus far, the influence of RT anneal- 
ing has been seen experimentally in the superconducting 
transition temperature, T c |Q, the Raman 2-magnon line 
(2-ML) and electronic continuum intensity 0, and var- 
ious structural || and optical Q] parameters. However, 
a detailed understanding of the kinetics involved has yet 
to be achieved. 

We have investigated the kinetics of several physical 
properties of underdoped YBa2Cu30g.4 on the same sam- 
ple during reproducible quench/anneal cycles, and find 
two characteristic behaviors. One is exhibited by the dc 
resistivity, T c , and the low-frequency (LF) optical spec- 
tral weight, all of which recover slowly over a few days 
following heat treatment. In sharp contrast, the opti- 
cal charge-transfer band (CTB), the B\ g Raman 2-ML, 
electronic continuum, and phonon intensity relax rapidly 
during the first hours of RT annealing. Since the dc re- 
sistivity and T c involve the Cu0 2 -plane carrier density 
directly, our results suggest that the rapid relaxation of 
the second group occurs through a process that is inde- 
pendent of chain-to-plane charge transfer. 

Monte Carlo simulations of oxygen ordering in the 
chains also show two relaxation rates: one fast, corre- 
sponding to short range oxygen order where the correla- 
tion length is of order the lattice constant, r « a; and one 
slow, where the correlation length is large compared to a. 
We conclude from these results that the fast relaxation 
processes depend on short range correlations, whereas 



long range oxygen order is important to the slower pro- 
cesses. 

As shown in the inset of Fig. [I], the underdoped 
YBa2Cu306.4 single crystal used for these measurements 
displays a sharp superconducting transition (AT C < 
1—1.5 K) at each quench/anneal cycle, despite T c changes 
from 18 K before the thermal treatment to 2.9 K immedi- 
ately following the quench. The reproducibility of these 
results for many cycles indicates that the total oxygen 
content of the sample remains constant. The homogene- 
ity of the oxygen distribution at the 1 fim level was con- 
firmed by micro Raman experiments. 

The optical reflectance at near-normal incidence in the 
energy range u = 0.03-2.25 eV was measured at succes- 
sive times following thermal treatment using a rapid-scan 
Michelson interferometer. A Kramers-Kronig transfor- 
mation || was applied to obtain the optical conductivity 
a(t,u>) from the measured reflectance data. Since varia- 
tions in the high frequency reflectance (3.8 < u < 4.3 eV) 
j| cause changes of less then 5% in the conductivity 
a(t,uj) for id < 2.5 eV, ellipsometric data measured af- 
ter complete relaxation were used in the energy range 
u> = 2.25-6 eV to make this transformation. The inset 
of Fig. H shows a(t, v) for YBa2Cu306.4 at two times 
t following the quench, to — 10 m i n and (fully re- 
laxed). Immediately after the quench, the largest changes 
in a(t, uj) are a decrease in the LF spectral weight (box 1 
in the inset), and an increase in the CTB spectral weight 
(box 2) J7J. The LF spectral weight (box 1) is comprised 
of a Drude part due to holes in the Cu02-planes plus a 
mid-IR peak associated with chain conductivity ||. In 
Fig. ^ we plot as a function of time the variation in ef- 
fective number of carriers (integrated spectral weight), 
AJVi )2 (f) = C f boxly2 [v{t,uj) - ff(tco,w)](iu, where the 
constant C is given in Eq. (3) of Ref. [7]. We find that 
AA^(i) takes ~ 100 min to recover to its initial value, 
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whereas complete recovery of AN±(t) requires several 
days. Moreover, the sum AN(t) = AJVi(t) + AN 2 (t) 
during the first ~ 100 min remains approximately con- 
stant. These results are significant because they suggest 
that the transfer of spectral weight from the CTB to LF 
band occurs quite rapidly relative to the time scale asso- 
ciated with chain-to-plane charge transfer. Notably, the 
monotonic increase in AAq(t) at long times may reflect 
different effective masses for holes in the chain and the 
plane layers. 

Raman scattering (RS) data from single crystal 
YBa 2 Cu306.4 was obtained in 20 second intervals at suc- 
cessive times following the quench using a CCD equipped 
Raman triple spectrometer. The 2-ML [||J9) and elec- 
tronic continuum ||,[lo| spectra were obtained at room 
temperature, while the 336 cm -1 phonon scattering data 
were taken at 5 K. The inset of Fig. || presents the RS 
spectrum measured at to — 6 min after the quench in the 
x'y' (mainly B\ g ) geometry. Fig. [3] shows the temporal 
dependence of the 2-ML and electronic background (0.1 
to 0.5 eV) intensities. An increase of both the 2-ML and 
low-frequency RS background intensities is observed af- 
ter quenching; however the fractional increase is larger 
for the background than for the 2-ML. Interestingly, the 
main change in the intensity for both features occurs dur- 
ing the initial 50-100 min, although complete recovery to 
pre-quench values takes several days. 

At least two characteristic relaxation rates may be 
extracted from Figs. 1-3 if an exponential temporal 
dependence of the form Ait) — Ait^) — [A(t ) — 
Aitoo)] exp(— f/r) is assumed. In order to estimate these 
times, we have represented the data as 
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where A{t) is one of the measurable quantities and A(to) 
is the first measured value of A just after the quenching. 

Figure ^a plots La{£) for the optical and RS quanti- 
ties, and T c . Two distinct relaxation behaviors are clearly 
seen. The kinetics of T c and AN±(i) demonstrate a sin- 
gle long relaxation time, estimated to be r w 250 min 
via a straight line least squares fit to the experimental 
points in the range t = — 200 min. In jllj], the kinetics 
of the dc resistivity and T c during room temperature an- 
nealing were shown to be of the same order. By contrast, 
the kinetics of 2-ML, background intensity, the integrated 
intensity of the B\ g phonon and o~(t, u>) in the CTB spec- 
tral region demonstrate two relaxation times: an initially 
short relaxation (r f» 50 min) over the first « 100 min, 
followed by a considerably slower relaxation. 

To gain insight into the experimental data, we carried 
out a Monte Carlo simulation of the oxygen distribu- 
tion kinetics in the CuO-chains of YBa2Cu30g+x- We 
used the asymmetric next-nearest-neighbor Ising model 
[|12|1 , which takes into account three different effective 



pair interactions between oxygen atoms (V! = 367 meV, 
V 2 = -0.348V! and V 3 = 0.159V!) @. This model has 
been used previously [ p^|JT^ ] for collecting snapshots of 
oxygen placement in the CuO-chains at each Monte Carlo 
step to study the kinetics of oxygen ordering, and also as 
a method for determining final oxygen structure Jl5| , ^6[ . 
We parameterize the results of these simulations in terms 
of an ordering parameter a, which measures how close the 
structure is to the ordered Ortho-I (Ol)-phase |l7| ]. a is 
defined to be zero when the structure is fully disordered 
and unity when the structure is in the ideal Ol-phase. 
We calculate a for different correlation distances r by 
comparing, at every oxygen position, the Monte Carlo 
oxygen arrangement on a square of side 2r to the perfect 
OI structure. We used up to 5000 Monte Carlo steps 
on a 128x128 lattice at room temperature, starting with 
a disordered structure. To compare the calculated and 
experimental data we defined function L a (t), for the or- 
dering parameter a that is similar to -La(^) (Eq. 1), where 
too is equal to 5000 Monte Carlo steps. Fig. ||b shows the 
time dependence of L a for different correlation distances 
r with oxygen content x — 0.4. It is seen that the depen- 
dence of L a on time for different r looks like the tempo- 
ral dependence of the various measured quantities. For 
long correlation distances L a demonstrates slow kinetics, 
whereas the time dependencies of L a for short distances 
(r = 1 - 2) demonstrate two rates: fast for the initial 
stage and much slower at longer times. The appearance 
of two relaxation rates at small r results from the fact 
that at the initial stage of relaxation only local ordering 
takes place. Consequently, many little twins (a few lat- 
tice constants in size) appear. This process is rather fast 
and the associated correlation lengths r are quite small. 
The slower kinetics correspond to motion of twin walls 
which results in the subduction of smaller twins by larger 
ones. This process is revealed at all correlation distances 
|l7| . Similar results hold for the Ortho-II ordering pa- 
rameter. As expected, calculations performed for higher 
values of the oxygen content show less difference between 
the kinetics at various correlation distances. 

The long relaxation times exhibited by the dc conduc- 
tivity fllJI , T c and AJVi(t) imply that the fast kinetic 
behavior displayed by the CTB and the B\ g RS features 
occurs before oxygen ordering in the chains causes sig- 
nificant hole transfer to the Cu02-planes. Since some 
rapidly annealing parameters, such as the 2-ML, are be- 
lieved to have their origins in the planes, we conclude 
that short range oxygen- ordering in the chains gives rise 
to short-length- scale changes in the planes without signif- 
icantly altering the average planar carrier density. 

Comparing Fig. |]a with Fig. ^b, one can conclude that 
the kinetics for AN 2 (t) and for the B\ g symmetry RS are 
determined by local oxygen ordering at nearest neighbor 
distances. As discussed in the CTB is most probably 
due to a local excitation of a hole in the Cu02-planes 
from Cu(2) d x 2_ y 2 to the four 0(2,3) p x , y states that 
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surround a first-neighbor Cu(2) site. The rapid transfer 
of the integrated spectral weight from AA^i) to AN±(t) 
suggests partial derealization of bound carriers. 

The 2-ML mainly probes short wavelength magnetic 
excitations from the vicinity of the magnetic Brillouin 
zone boundary k = (7r/a,0), which occur via photon- 
driven superexchange between two nearest-neighboring 
Cu(2) spins in the antiferromagnetic surroundings. Thus, 
the CTB and 2-ML both require only short range anti- 
ferromagnetic order which is affected by changes in the lo- 
cal carrier distribution but is independent of the average 
carrier number. The rather fast kinetics of the Raman 
Big continuum implies that it is also a local excitation, 
perhaps of magnetic origin, and that it is not related di- 
rectly to the average carrier number. Phonon interaction 
with the continuum might explain the observed phonon 
line intensity kinetics through a Fano-type effect, but a 
detailed model needs to be formulated to understand this. 
By contrast, the resistivity, T c , and the LF absorption ki- 
netics are determined by longer distance oxygen ordering 
connected with the formation of ordered oxygen chains 
several lattice constants long. The LF band has been 
shown to be polarized along the chain direction ; thus 
chain fragments of a minimum length are required for 
this band to occur. 

In summary, we have measured the kinetic behavior 
of the magnetic susceptibility, the B\ g symmetry Raman 
scattering (2-ML, electronic continuum and the 336 cm -1 
phonon), and the optical conductivity (LF and CTB) 
of YBa2Cu3C>6.4 single crystals during room tempera- 
ture annealing. The kinetics of the Raman features and 
the optical conductivity in the CTB region demonstrate 
fast relaxation (ry ~ 50 min). The relaxation time for 
the resistivity, the LF absorption, and T c is much larger 
(r s > 220 min), implying that there is a nearly constant 
average carrier number for times less than 77. Compar- 
ing the measured kinetics with Monte Carlo model simu- 
lations we suggest that the short-length-scale excitations 
arc responsible for the fast processes, such as B\ g symme- 
try RS and the CTB in the optical conductivity. Longer 
correlations (r > 2) are responsible for changing the re- 
sistivity, LF optical conductivity and hence the effective 
carrier concentration. 
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FIG. 1. Temporal dependencies of T c after the quench. 
Inset: Dependence of magnetic susceptibility \ on tempera- 
ture measured at different times after the quench. 

FIG. 2. Temporal dependencies after the quench of the 
variation in the optical conductivity integrated over the LF 
(open circles) and CTB (filled circles) regions. Inset: Optical 
conductivity a(t, ui) before the heat treatment and 10 min 
after the quench. Boxes 1 and 2 show the integration regions. 

FIG. 3. Temporal dependencies of the intensity of the 
2-magnon Raman scattering line maximum (open circles) 
and the background intensity (filled circles) after the quench. 
The vertical axis scales are anchored to the time axis at 
their pre-quench values (i.e., 7.7 for background and 38.5 for 
2-ML). Inset: Raman spectrum measured in x'y' geometry 
at t = 6 min after the quench. Filled region corresponds to 
broad structureless electronic background. 

FIG. 4. (a) - Measured values La(^) [Eq. (1)] for parame- 
ters in Figs. 1-3. Straight lines are the results of least squares 
fitting of experimental data at the initial stage of aging. 
Crosses - T c , r ~ 340 min; squares - AiVi(t), r « 220 min; 
filled circles - AN2(t), r ~ 59 min; triangles - background, 
r « 51 min; open circles - 2-magnon line, r « 36 min. 

(b) - Time dependencies of L a for calculated pa- 
rameter a in YBa2Cu306.4 crystals for different correlation 
distances 1 < r < 6 in units of \/2a. 
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